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ABSTRACT 



We consider the motion of a charge in a large amphtude electrostatic wave 
with a triangular wave form relevant to an oscillating model of a pulsar magneto- 
sphere. The (one-dimensional) orbit of a particle in such a wave is found exactly 
in terms of Weierstrass functions. 

The result is used to discuss linear acceleration emission (at both low and high 
frequencies) in an oscillating model for pulsars. An explicit express ion for the 



20091), and used 



emissivity is derived in an accompanying paper ( iMelrose fc Lud 
here to derive an expression for the absorption coefficient at low frequencies. We 
show that absorption can be negative, corresponding to maser emission. For the 
large amplitude required to trigger pair creation in an oscillating model, the rate 
of the maser growth is too small to be effective. Effective growth requires smaller 
amplitudes, such that the maximum Lorentz factor gained by acceleration in the 
wave is < 10. 



Subject headings: plasmas — pulsar: general — radiation mechanism: nonthermal 
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Introduction 



An oscillatory model for a pulsar magnetosphere introduces an emission mechanism 
that depends explicitly on the oscillating electric field: line ar acceleration emissio n (LAE). 



20091), referred 



Our purpose in this paper, and in an accompanying paper (iMelrose &: Lud 
to hereinafter as paper 2, is to develop the theory for LAE in a large amplitude electrostatic 
wave (LAEW), and discuss its viability as a pulsar emission mechanism. 

In polar-cap models for pulsars, the inductive electric field associated with the 
rotation of the magnetized neutron star cannot be completely screened everywhere, 
and the residual field in 'gaps' leads to acceleration of particles along the field lines 
to high enough energies to trigger a pair cascade. The polar-cap region is populated 
by the resulting secondary pairs, which escape through the light cylinder to form the 
pulsar wind. In most models, the electric field, the acceleration, the pair creation, and 
the screening of the inductive field are assumed to occur in a time-stationary way in 
a frame corotating with the star. An alternative view is that these processes occur in 



an oscillatory way fISturrock 



1971 



Al'Ber. Krotova fc Eidman 



19751 ). Recent numerical 



models confirm that the system is unstable to temporal perturbatio ns that develop into 



large amplitude electrostatic oscillation s, both for ordinary pulsars ( 



Melrose. Levinson. Judge, fc Luo 



2007a 



Levinson et al 



2005 



20051 ) and for magnet ars (IBeloborodov. fc Thompson 



b|). In these nurnerical models the oscillations are purely electrostatic and purely 



temporal. 



oscillations, and 



Beloborodovl (120081) p r opose d Alfven waves as an alternative to the electrostatic 



Luo fc Melrosd (120081 ) proposed an analytic model that generalized the 



purely temporal oscillations to large amplitude propagating waves. The analytic model 
includes both purely temporal an d purely spatia l oscillations as limiting cases; purely 



spatial oscillations were found by 



Shibatal (119971 ). More generally, the analytic model 



describes an outward-propagating, large-amplitude electrostatic wave; the phase velocity. 
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cl3y, is a free parameter with purely temporal oscillations corresponding to /?y oo, and 
purely spatial oscillations to /5y = 0. A bulk outward motion occurs at the counterpart of 
the group velocity, which corresponds to a speed c/ [3v in the superluminal case, [3y > 1. 
For > 1 one can make a Lorentz transformation to the frame moving at (3* = —l/(3v, 
referred to here as the primed frame, in which the oscillations are purely temporal, and it is 
convenient to adopt this frame for detailed calculations. We envisage a model for a pulsar 
magnetosphere in which the pair creation occurs in sporadic, isolated events association 
with the development of LAEWs. It is possible that such a model, when averaged over 
the effects of the LAEWs, mimics a conventional time-stationary model, but it is more 
likely that new features introduced by the LAEWs result in important differences from 
conventional time-stationary models, with the emission of LAE being one such new feature. 

Pair creation is included in the os cillating model, but has not been treated in the 



same detail as in the stationary model (jZhang fc Harding 



Arendt fc Eilek 



2000 



Hibschman &: Arons 



2001 



2OO2I ). The rate of pair creation incre ases with the amplitu de of the LAEW 



2OO5I ). We assume 



providing a mass loading that limits the amplitude (iLevinson et al. 
that the pair creation produces particles with a spread in Lorentz factors that is small 
compared with the maximum Lorentz factor, 7max, achieved in the LAEW, and we neglect 
this spread for the 'background' particles that are described in terms of their (oscillating) 
velocity, Lorentz factor and number density. The results derived here are not sensitive to 



the frequency o : 



( iLevinson et al. 



the LAEW, which we estimate to be ~ 10*^-10^ s ^, based on ujp/'y 



1/2 



2OO5I), with 7n 



10^-10^ and the plasma frequency determined by the 



Goldreich- Julian number density times a multiplicity M: Up ^ ^/2^l^^lM ~ 5 x 10^^ s ^ for 
M = 100, a cyclotron frequency, fig for B = 10^ T, and a rotation frequency, Q = 2tt/P 
with P = 0.1 s. 



In the present paper, we discuss the motion of a charge in a LAEW, and the associated 
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emission of radiation. A first integral of the equation of motion is obtained for a particle 
in a LAEW with an arbitrary wave form. A second integral requ ires a specific choice of 



wave form. The solution for a sinusoidal wave form was found by 



Rowd fll992a 



b|) in terms 



of elliptic functions. Here we consider a sawtooth or triangul ar wave form, w . 



lich is an 



excellent approximation to the actual wave for m for a LAEW 



case where the particles are highly relativistic (iLuo fc Melrose 



Luo 



fc Melrose 



20081) in the 



20081 ). The solution for the 



orbit is found in terms of Weierstrass elliptic integrals. The exact solution is used as the 
basis for a treatment of LAE in a LAEW, which is developed in terms of Airy integrals in 
paper 2, and used here to discuss the possibility of maser LAE in a LAEW as a possible 
pulsar radio emission mechanism. 

The equation of motion for a relativistic particle is solved in two superficially different 
ways. In Appendix [Xj a covariant form for the wave equation is used to derive a solution 
in terms of (Lorentz) invariants; the solution may be written in any inertial frame by 
expressing the invariants in terms of quantities in that frame. In Sec. |2] the solution is 
obtained in the primed frame, defined to be the frame in which the oscillations in the 
LAEW are purely temporal. (There is only a notational difference between the important 
results in Appendix |A] and Sec. [2l) The exact orbit for a triangular wave form is found in 
Sec. [HI The application to LAE is discussed in Sec. HI concentrating on the possibility of 
maser LAE at low frequencies. The results are discussed in Sec. |5l and our conclusions are 
summarized in Sec. [61 



2. Motion in a LAEW: primed frame 

In this section we find the orbit of a charge in a LAEW in the primed frame. First, we 
describe the transformation between the inertial frames, and write down the forms of the 
invariants introduced in Appendix [X] in the primed frame. The motion is assumed to be 



restricted to one dimension (the 3-axis), along the direction of the superstrong magnetic 
field in a pulsar magnetosphere. 



2.1. Description of the LAEW 

The LAEW is described in terms of its amplitude, E{x) = EqT[x)i is a function of 
phase, X, with the maximum electric field, Eq, written in terms of a frequency, uje, and 
with the wave form described by the function T(x), with maxima and minima ±1: 
77 — > e everywhere 

E[x) = e MX , t^fi = , 1 

q mc 

where the sign, e = q/Iq], is opposite for electrons and positrons. For an arbitrary wave 
form, the only specific assumption made is that T{x) = T{x + 27r) is periodic. A LAEW, 
even with a non-sinusoidal wave profile, is described in terms of its frequency, fl, and its 
wave vector, K, along the 3-axis. The LAEW is assumed superluminal, implying that its 
phase velocity Vl/K = j3vc, satisfies /3y > 1. 



2.2. Lorentz transformation 

The Lorentz transformation to the primed frame corresponds to a boost with velocity 
(in units of c) (3* — —l/Pv and Lorentz factor 7* — PviPv ~ The relation between 

time and position in the two frames is 

ct' ^^*{ct- p*z), z' ^f{z- (5*ct). (2) 

In the primed frame, the frequency, Vt\ of the oscillation and the 4- velocity, u^' = 7'[1, 0, 0, j3'] 
of the test charge are given by 
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where the wave vector in the primed frame is zero by construction, K' = 0. The frequency 
and angle of emission of LAE are related in the two frames by 

. = ,V(l+;3-cos.O. cos» = ^^^. (4) 

Relevant invariants introduced in Appendix |A] have the following values in the primed 
frame, implied by K^' = (^]', 0, 0, 0), JsT^' = (0, 0, 0, n'): 

n' = {K^f/\ x = n't', Ku{x) = n'i, Kj,u{x) = -n'u\ (5) 

where we assume = — |Kp > and write u' = 7'/?'. The electric field, E{x), is an 
invariant, and has the same form in all frames. 



2.3. Equation of motion 

The equation of motion may be written in two equivalent forms: 

Integrating the second of gives 

u'{x)=u', + {u;E/n')F{x), (7) 

where u'q = 7q/5q contains the initial conditions, and with F{x) given by flASI) . with oq = 
assumed in (lAlOp . The solution ([7j) implies a Lorentz factor 

7'(x) = {l + K + (^i./Wx)f}'^ (8) 
and a velocity (in units of c) 

3'M = < + Mn')F{x) 

{i + [u', + {ujE/n')F{xW}y^- 



- 8 - 



2.4. Background and test charges 

In the solution ([7]), the particle has an arbitrary initial velocity, /3qC, at x = 0. It is 
convenient to separate the particles into two classes: background particles and test charges. 
The background particles are an idealized class of electrons and positrons that move in 
opposite directions in phase with the LAEW, with /3q assumed to be identically zero. In the 
idealized case where there is no intrinsic velocity spread in the background particles, such 
particles are instantaneously at rest twice per period, and reach their maximum Lorentz 
factor twice per period. Let the initial phase, x = 0, be chosen such that the background 
particles have /3q = in ([7]), so that they are instantaneously at rest ai x = ^'^ for any 
integer n. The background particles have their maximum Lorentz factors, 7max5 ^^e 
phases x = 7r/2 + n7r. A test charge is then any particle that has (3'^ ^ at x = 0- In the 



analytic model (ILuo fc Melrose 



20081 ) all particles are in effect test charges. The concept of 
a background particle is useful when considering situations in which the velocity spreads of 
the electrons or of the positrons are unimportant, and are approximated by zero. 

The maximum 4-velocity for the background particles, w^ax — 7max/^maxi determined 
by dZ]), which becomes 

The Lorentz factor corresponding to ( ITOl) is 

7'{x) = {l + {uJE/nr[F{x)?V^'. (11) 



2.5. Orbit 

A second integration gives the orbit: 



z' = z', + ^Z\x\ Z\x) = £dx'f3'{x'). (12) 



- 9 - 



The explicit form for Z'{x) obtained by inserting the solution iQ into f|T2l) cannot be 
evaluated in closed form in general. It is evaluated for a triangular wave form in Sec. O 

The motion implied by (fT2l) for a background particle, Mq = 0, is an oscillation around 
a fixed point; the particle returns to the same position at the same phase in each subsequent 
oscillation, and its mean velocity (and 4- velocity), averaged over a period of the LAEW, is 
zero. A test charge, u'q ^ 0, has a net drift motion, in the sense that its velocity (and its 
4-velocity) averaged over a period is nonzero. One may regard a test charge as oscillating 
about a drifting point. For 7q < 7^ax) ^ ^^st charge is instantaneously at rest twice per 
period; for 7q > 7^ax5 P'ix) never changes sign. 



3. Triangular wave form 

For oje/^' ^ 1 the electric field is strong enough to cause all particles to be highly 
relativistic for all phases of the LAEW with the exception of a brief phase where the 
direction of their velocity reverse s. The LAEW is then well approximated by the triangular 



wave form (ILuo fc Melrose 



20081 ). illustrated in Fig. [H An interpretation of the triangular 
wave form follows by noting that the current density, J, associated with relativistically 
counterstreaming electrons and positrons is approximately unchanged by the acceleration 
(by E parallel to J), and that the solution of dE/dt = — J/^o implies that E varies linearly 
with t] J reverses sign abruptly in the brief phases where the particles are nonrelativistic, 
and E reverses sign out of phase with J. Here we solve explicitly for the orbit for such a 
triangular wave form. 
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3.1. Form for T{x) and F{x) 



Our model for a triangular wave form is 



2x/7i-l, 0<x<vr, 
T{X) = { (13) 
3-2x/vr, 7r<x<27r. 



The function -F(x), defined by (lASp . is 



i^(x) = - <^ (14) 
^ ' -(X - 37r/4)2 + 7rV4, tt < x < 27r. 

A background particle has extrema in its velocity at x = 7r/2 + nn, so that ffTOl) implies 

KaxI = T^UJE/^n', (15) 

with I Umax I ~ 7max ^ 1 in the case of most interest here. 



3.2. Exact solution for the trajectory 



The integral over x' in ([12]) with may be carried out exactly for the triangular wave 
form, in terms of elliptic inte grals, specifically, the Weierstrass p{x) and ({x) functions 



(lAbramowitz fc Stegun 



19651 ). The integral is written as a sum over integrals of the form 
(X - a)' 



dx 



P 



(16) 



where xi, X25 and h are constants. By making the substitution y = {x~ '^Y ~ 2p/3, the 
integral ( fT6l) reduces to the form 



± 



y2 



(17) 



4(v'-3b)/3, Q, = - 


-8p 


[p^ + 


(Abramowitz & Stet 


run 


1965) 



A/4y'3 - g2y' - 5(3' 
96) /27, which is a standard form for the Weierstrass 



C(a;;5'2,5'3) = - / dxp{x; 02,93), 
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denoted ({x) here, with p{x; g2, gs), denoted p{x) here, defined imphcitly hj y = p{x) and 



X 



dy' 



(19) 



Two sets of choices, ai,pi with i = 1,2 are needed to find the full solution, and these 
define 

x^{x) = p-\{x~a,Y~2p,/3). (20) 



In terms of the functions 



Wiix) = PiXiix) + K{xi{x))i 



with (72i, g^i evaluated for p = pi, the solution is 



z'ix) 



-w^{0) + wi{x), 0<x<vr/2, 

Wi{tc/2) - wiix)- 7r/2<x<vr, 

^2(71-) - W2(x), 7r<x<37r/2, 

-«;2(37r/2)+W2(x), 37r/2 < x < 27r, 



(21) 



(22) 



with ai = 7r/2, pi = n^/A - 7t%Po/u'^^^, = 37r/2, p2 = 7rV4 + n%p'Ju'^^^, b = n^/u[ 



2 /„,/2 

max' 



3.3. Examples 

Examples of the exact solution are shown in Fig. O A feature of the solutions is that the 
trajectory is well approximated by straight lines, with one of z'^ct' constant, corresponding 
to Z'{x) ~ ill the first half phase, with Z' oscillating between approximately ±7r. The 
linear approximation is accurate except for a small range of phase, ~ l/7max5 where the 
particles become nonrelativistic and reverse their direction of propagation. As is evident 
from Fig. [21 the trajectory of a background particle (/^q = 0) involves an oscillation about a 
fixed point in the primed frame. A test particle with ImqI < iM^axI; reverses its direction of 
motion periodically in the LAEW, but for ImqI > I'^maxl' ^^e direction of motion of the test 
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0.5 



-0.5 




2n 



3 7T 
1 



4 7T 



5 7T 



6 7T 



i^'ig. 1. — The trian gular wave form is illustrated for a LAEW; in the exact theory 



Luo fc Melrose 



20081 ). the sharply curved parts correspond to a range of phase ~ 1/7^ 



where the background particles are briefly nonrelativistic. 




Fig. 2. — The trajectory of a background particle (solid curve for an electron, dashed curve 
for a positron) in a LAEW for a modest relativistic value, "Umax = 10 (7max = 10.05). For a 
more realistic value, 7max ~ ^O^; ^^e orbit is very closely approximated by a triangular form, 
consisting of a sequence of oppositely directed light lines. 
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particle does not change, as its 4-velocity oscillates in the LAEW. For \u'q\ ^ l^^maxl the 
effect of the LAEW on the motion of the test charge may be treated as a perturbation to 
the motion of a background particle. 

For 7^j^^ ^ 1, the orbit of any particle (except for |7max ~7ol ^ 1) ^^^^ approximated 
by a sequence of light lines, such that z' varies as ±ct' in sequential sections, as shown on 
the left in Fig. [31 For the large value, 7^^^^ ~ 10^5 expected in an oscillating pulsar model, 
the fraction of the phase. Ax ~ l/7max5 where the particles are nonrelativistic is tiny, and 
the approximation in terms of a sequence of light lines is very accurate. The distance 
between the turn-around points for a background particle is very nearly the light distance, 
nc/Q', corresponding to the half period in which the particle travels in the same direction. 
The solution in terms of Weierstrass functions shows how the orbit evolves from one section 
to another; to illustrate this, an example where the motion is only mildly relativistic is 
shown on the right in Fig. [3l 

4. Particle current associated with LAE 

The motivation for the investigation reported is the treatment of LAE due to motion 
in a LAEW. In Appendix |X] we use the covariant formalism to evaluate the 4-current 
associated with LAE, and in this section we use the results derived in Appendix [A] to write 
down the 3-current associated with LAE in the primed frame. 

4.1. LAE in the primed frame 



The emission of LAE can be calculated in terms of the current associated with the 
motion of the ch arge in the LAEW. For periodic motion, the current may be expanded in 
a Fourier series (lRowelll995l ). giving a sum over a harmonic number s. The expansion is 
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written down using a covariant formalism in Appendix \M In the primed frame, the current 
along the direction defined by the LAEW follows from (IA22p : 



Jsiiu', k') = -Bs{iu', k') 2n5{iu' - sn'), (23) 

with Bs{uj',k.') following from flA20p . It follows that in the primed frame, emission at the 
sth harmonic is at the frequency, u' = sQ'. That is, the emission is at harmonics of the 
frequency of the LAEW in the primed frame. 

The explicit form for the current obtained by inserting (lA20p into fl23p involves an 
integral over a phase factor, with phase (f){x) say. After replacing the sum over s by an 
integral and performing the integral over the 5-function in fl23l) . the phase becomes 

<P{x) = ^[x-Z'{x)cose'], (24) 

in the primed frame. 

On inserting the exact solution (l22l) for the orbit into f l24p . it is obvious that some 
simplifying procedure is needed before one can evaluate the integral over the phase factor. 
The analogy between LAE and synchrotron emission, discussed in detail in paper 2, suggests 
that an Airy-integral approximation is appropriate. The Airy integral applies when the 
phase is approximated by an expansion about a specific phase in which the quadratic term 
vanishes and only the linear and cubic terms are retained. On making a Taylor expansion 
of Z'{x) about an arbitrary phase, Xo say, the quadratic term vanishes when the second 
derivative of Z'{x) is zero at x = Xo- This corresponds to an extremum of /5'(x) and hence 
of 7'(x). The only relevant extrema are maxima of 7'(x)) which occur at = + nir. 
Assuming the radiating particle to be highly relativistic near this phase, one may make the 
approximation l3'{x) ~ =^[1 ~ V27'^(x)]- This approximation leads to a highly relativistic 
counterpart of fl22|) . For Uq = the integral can be performed exactly, giving 

l + 2(x-7r/2)/7r 



2(x-7r/2)/7r 



(25) 
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for < X < TT. Comparison with fl22|) shows that fl25l) is an excellent approximation except 
near x = 0, vr where the particle is briefly nonrelativistic. The Airy-integral approximation 
is obtained by expanding fl2^ around Xo = ^/2- The same result is obtained by expanding 
7'(x) around it/2. The latter procedure is readily generalized to u'q ^ 0, for which the 
maxima in 7'(x) also occur at xo ~ 7r/2 + nvr and are ~ 7^ax i ""o- The Airy-integral 
form for the phase fl24l) then follows from 



z'ix) =x-7^ 



_ ^^ 8 7max / _ VTAjS' 

2/ ^3^^V^^ 2 



(26) 



27i L 

The evaluation of the integral is discussed in detail in paper 2, and here we quote the result 
and discuss some of its implications. 



4.2. Emission and absorption of LAE 

The formal treatment of LAE involves both conceptual and mathematical complications, 
discussed in paper 2. To proceed with the discussion here, we quote the results derived in 
paper 2 for the emissivity (power per unit frequency and per unit solid angle) in the primed 
frame. In terms of the characteristic frequency = 2f2'7^, one has 

V±iu;,e)- L±Ky,{^±), (27) 

where Ki/3{^) is a modified Bessel function, with = ^c±{^ + 7^6*'^)^^^, ^c± = 
{n / ^/2){lu' /uj'^j_){'y'^/^'^^^y^'^ , and where the ± signs correspond to emission in the forward 
(defined for u'q > 0) and backward direction. The characteristic properties of LAE implied 
by (|27|) are discussed in paper 2. 

Let g'{u'Q)du'Q be the number density of one species of particle, electrons say, with 
initial (at phase x = 0) 4- velocity in the range u'q to u'q + du'^. The volume emissivity 
is J duQg'{u'f))ri'^{uj' ,6'). The absorption coefficient can be related to the emissivity (!27|) 
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by an argument using detailed balance (lTwisall958l : IWild. Smerd fc Weisall963l ). In the 
one-dimensional case of relevance here this gives 



TW.e') = / du',v'A^',e') cos^'^^^^'o 



duQ 



(28) 



with cosO' ~ ±1 such that the signs of cos^' and Uq are the same. Maser emission 
corresponds to negative absorption, requiring r'^{uj',6') < 0. We comment below on the 
possibility of maser emission at low frequencies. 

As discussed in paper 2, a mathematical inconsistency arises when one attempts to 
consider the low- frequency limit of fl28l) with fl27|l . One way of avoiding this difficulty is to 
consider the average of fl28l) over solid angle, denoted by an overline: 



r±(^') 



/ duQr]^{u:) 



muj ■ 



duQ ' 



(29) 



where the choice of the sign cos^' = ±1 is implicit. The average of the emissivity is 
evaluated in paper 2, and at low frequencies for a triangular wave form it reduces to 



2q^n' fnco' 



4/3 / \ 2/3 



Ai2(0) 



(30) 



with Ai(0) = l/32/3r(2/3) = 0.355. We use ^ rather than ([28]) in the discussion below. 

The emissivity f l27l) is in the primed frame, and it is straightforward to transform it 
to the pulsar frame, in which the LAEW is an outward propagating wave. The Lorentz 
transformation implies 



UJ 



7*(l±/5*)cu', e = -i*{lTl3*)e' 



(31) 



However, it is usually more convenient to discuss the emission and absorption in the primed 
frame before transforming to the pulsar frame. 
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5. Discussion 

In this section we discuss four aspects of tlie foregoing results in further detail: 
systematic acceleration, the displacement of individual charges in a LAEW, drift motion 
and its relevance to LAE in a LAEW, and maser LAE in a LAEW. 



5.1. Systematic acceleration 



In solving the equation of motion, we find that a LAEW with an asymmetric wave 
form, specifically one with 7^ in flAlOp . leads to systematic acceleration of all particles. 
According to (jTj) with (lAlOp . the 4- velocity increasing by 27rao each wave period. This 
corresponds to a systematic acceleration by a static electric fi eld E = anfl'mc/q. N o uniform 



20081 ). and we 



electric field is included in the analytic model for a LAEW (ILuo fc Melrose 
assume oq = in the foregoing. 

It is interesting to speculate on the possible significance of assuming 7^ 0. One could 
interpret the E associated with 7^ as one way of modeling an unscreened component 
of the pulsar's inductive field. The systematic acceleration transfers electromagnetic energy 
into particles. In a conventional polar-cap model, such acceleration occurs in a gap, 
and s creening of the inductive fie ld results from a net charge density in a pair formation 



front (IHarding fc Muslimov 



19981 ). One motivation for an oscillating model is that such 



screening is unstable to tempora l perturbations, resulting in large-amplitude electric 



oscillations flLevinson et al. 



20051 ). The build up of a LAEW in an oscillating model relies 
on acceleration by an incompletely screened electric field. As in a stationary model, the 
resulting pair creation should lead to screening of the electric field, and whereas this occurs 
locally in a stationary model, it corresponds to a systematic reduction in E, and hence of 
Oo, in the oscillating model. This effect needs to be included in a detailed theory of the 
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instability leading to the LAEW, but we do not attempt a quantitative treatment of this 
here. 



A particle in a LAEW oscillates about a center that is drifting (except for a background 
particle in the primed frame). The displacement about this center is by ivrc/fi' in the 
primed frame. A particle is constrained to move along a magnetic field line that is curved, 
and the neglect of this curvature is valid only if this distance is small compared with the 
radius of curvature of the field line. For a LAEW with a frequency ~ 10^ s^^, ±7rc/fi' is less 
than a meter, and this condition is well satisfied. 

A test particle has a drift velocity in the primed frame. This velocity may be estimated 
from the distance, I\z' say, its position advances in each period of the LAEW: 



with given by ([9]). The mean drift velocity is i7'A2;'/27r, which depends on Mq, and 

Wmax, and is zero for Mq = 0. The 4- velocity constructed from the 3- velocity Vt' Az' /2tt is 
rather cumbersome. Alternatively, one might consider a drift 4-velocity found by dividing 
the displacement (15^ by the proper time elapsed in a wave period, Ar say. One has 



with 7'(x) given by ([8]). The ratio Az'/At is some measure of the mean 4-velocity, but 
describes a different quantity. We note that neither the integral fl32|) nor the integral fl33l) 
appear in the detailed theory for LAE developed in paper 2. 



5.2. Displacement of charges in a LAEW 




(32) 




(33) 



- 19 - 



5.3. Comparison with iRowd (1l995l ) 



The solution found here for the orbit of a part i cle in a LAEW with a triangular 



Rowd (I1995I ) for a sinusoidal wave form. 



ts for 



AE. An important difference 



wave form is similar to the solution found by 
However, there are important differences in the resul 
is the dependence of the emission formula found by 
we comment on this below. First, however, we comment on similarities in the results, to 
eliminate them as possible reasons for the differences in the treatment of LAE. 



Rowd (I1995I ) on a drift velocity, and 



As in the present paper, iRowd (Il995l ) found that when the particles are highly 
relativistic over most of the phase, their orbit is well approximated by a sequence of light 
lines, with the directi o n of m otion reversing at the phase where the particles are briefly 



nonrelativistic. iRowd (jl995l ) also found that the emission of LAE is dominated by the 
phase where the particles have their maximum Lorentz factor. From one perspective this 
is a surprising result, which is clearly insensitive to the assumed wave form. Based on the 
generalized Larmor formula, one expects the emission to be maximum when the acceleration 
is a maximum, that is, at the phase where the electric field is maximum. In contrast, it 
is found that the emission of LAE maximizes around the phase where the acceleration is 
zero. Clearl y, the i mpor tant differences in the treatment of LAE between the present paper 



and that of 



Rowd (I1995I ) are not due to the assumed wave form, to the detailed form of the 



orbit or to the phase that determines the properties of the emission. 



5.4. Absence of Doppler shift in LAE 



In the theory of LAE developed by iRowd (Il995l ). the drift velocity plays an 
important role, in the sense that the frequency of LAE is at ha rmoni c s of t he Doppler- 
shifted frequency of the LAEW. The emission in the theory of iRowd (119951 ) occurs at 
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[uj — kzcPo) — s{Q ~ KcI3d) = 0, ot kuo — sKud in invariant form, where is the drift 
4-velocity, 7d, I3dc are the associated Lorentz factor and 3-velocity, and s is the harmonic 
number. In contrast we find no exphcit dependence on a drift velocity in our treatment of 
LAE. 

ogous theory for emiss i on due to 



1971 



Arons 



mm. 



A drift velocity does appear naturally in the ana 
motion in a large amplitude transverse wave (LATW) (iGunn fc Ostriker 
In a LATW the oscillatory motion is in the transverse plane and there is a drift motion 
along the direction of wave propagation. The frequency of emission by such a particle is 
at harmonics of the Doppler-shifted wave frequency, with the Doppler shift corresponding 
to the drift velocity. However, the analogy between the transverse and longitudinal cases 
is misleading when considering the emission of radiation. Semi-quantitatively, one may 
regard the motion in a LATW as consisting of two effectively independent motions: that in 
the transverse plane causing the radiation, and that along the axis providing the Doppler 
shift. In a LAEW both the drift and the oscillatory motion are along the axis, and the two 
motions cannot be separated in this way. 

No Doppler shift associated with the drift motion of a test particle appears explicitly 
in the phase (12^ . As shown in Appendix |X1 explicit inclusion of the drift motion involves 
extracting a term, aix from Z'{x)^ such t hat Z'jy) = Z'{x) " C'lX is periodic. This does 



introduce a Doppler shift, as assumed by (IRowe 



19951 ) ■ In the continuum approximation, 



when s is replaced by a continuous variable, the sum over s is replaced by an integral that 
is performed over the resonance condition, s is replaced by the Doppler shifted frequency 
divided by the frequency of the LAEW. The phase is given by irrespective of whether 
the Doppler shift is included or not. If the Doppler shift is ignored, one has s = u'/Q' 
and fl24p follows by the argument given; if the Doppler shift is included, the term aix is 
subtracted from Z'{x)j added to s, and restored to Z'{x) to give Z'{x) after integration of s. 
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When the integral of the phase factor over x is evaluated in terms of an Airy integral, as in 
paper 2, no separation of the term aix associated with the drift is necessary or appropriate. 
It follows that no Doppler shift appears explicitly in LAE. 

A physical explanation for the absence of an explicit Doppler shift in LAE (in the 
primed frame) is as follows. The emission of LAE is concentrated around the phase where 
7'(x) reaches a maximum, with d'-f'{x)/dx = 0- This is the phase where the acceleration 
reverses sign and is instantaneously zero; all particles have their maximum Lorentz factor 
at the same phase. The theory of LAE depends only on the maximum value of 7'(x) ^^^^ 
the way that it varies with x around this maximum. The emission does depend on the 
initial 4- velocity, Uq, which also determines the drift motion, but the two dependences on 
u'q are quite different. In the emission of LAE, u'q affects the shape of the variation of 7'(x) 
about its maximum. In contrast, the drift velocity is determined by the integral (|32l) for the 
displacement in a period of the LAEW, and this is dominated by the asymmetry between 
the forward and backward paths, as shown in Fig. [3l in this affects the turn-around 

points, where 7'(x) passes through its minimum of unity. The dependence of LAE on Mq 
cannot be described in terms of the drift velocity. 



5.5. Maser LAE at low frequencies 

Maser LAE is known to be possible for small amplitude oscillations, when the frequ ency 



of th e maser is determined by a characteristic frequency in the radio range (IMelrose 



19781 ). Here the characteristic frequency, u'^^ ~ 2Q'y^^^, is at a much higher frequency, 
which is estimated in paper 2 to be in the hard X-ray range for the largest-amplitude 
LAEWs. For maser LAE in a LAEW to be relevant as a pulsar emission radio emission 
mechanism, the maser needs to operate at low frequencies, u ^ uJc±. It follows from 
fl28|) that negative absorption requires two conditions be satisfied. First, the distribution 
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function must satisfy dg{uQ)/duQ > 0, over at least some range of Uq. This condition is 



plausibly satisfied 



( iZhang &: Harding 



'or pairs generated in a cascade . Numerical models fo r pair creation 



2000 



Hibschman fc Arons 



2001 



Arendt fc Eilek 



20021) suggest that the 
distribution of particles increases with energy with a peak between 7q of a few to ~ 10^. 
**These calculations are for a stationary model, and we assume that the results are similar 
for the oscillating model, such that the intrinsic spread in initial Lorentz factors is ~ 10^. 
This spread is small compared with 7max ~ 10^-10^. Due to this spread, it is plausible 
that the initial distribution function is an increasing function of Uq with a peak in the 
range 1-10^, implying dg{u'f^ / du'^ > in fl25]) .** Second, by partially integrating fl2S]) . one 
finds that dri'^{uj' , 9') / du'^ < is also required for negative absorption to be possible. As 
already remarked, applying this condition to the emissivity itself encounters inconsistencies, 
and we avoid these by considering the angle-averaged form fl29l) . Using (150]) . one finds 
ri'j^{(jj') oc l/7±, and this requirement is well satisfied. 

The growth rate for LAE follows by partially integrating fl29|) and using (|30|) : 



r±(^') 



2(27r)^c 

^ / duQ g [uq 



muo' 



7± 



(34) 



where \d'^y'J^/ du'Q\ = 1 is used. Assuming u'q is negligible in comparison with 7^ax5 finds 



7/3 /^'\5/3 



21/3 



UJ' 



In 



8/3 



(35) 



The theory is valid only for u' ^ Q' and 7max ^ 1, so that both factors on the right hand 
side are very small. It follows that although the condition for growth at low frequency 
are well satisfied, the growth rate is very small. To explain pulsar radio emission in terms 
of maser LAE in a LAEW one would need LAEWs with much smaller amplitudes, with 



7; 



10. 
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6. Conclusions 

In this paper, we solve the equation of motion for a charged particle in a large 
amplitude electrostatic wave (LAEW), both in an arbitrary frame (Appendix and in the 
(primed) frame in which the oscillations are purely temporal (Sec. The case of interest 
for pulsars is when the electrons and positrons are accelerated to highly relativistic energies 
in a small fraction of a period of the wave, and in this case the wave form of the LAEW 
is well approximated by the triangular form Fig. [TJ The orbit of a particle is described in 
terms of its displacement, z, along the axis defined by the LAEW, as a function of phase 
X- We choose the zero of the phase such that the maxima and minima of the wave form for 
the electric amplitude at % = nvr, with n an integer. 

We distinguish between background particles and test particles. Background particles 
are assumed to come to rest at the maxima and minima in the wave form, and reach a 
common maximum Lorentz factor, 7^ax intermediate phases, x = + uti. A 

test particle is any particle with a nonzero 4-velocity (mq = 7o/5o, 7o = (1 — /5o^)~^^^) at 
the initial phase, x = 0. The Lorentz factor of a test charge oscillates between maxima 
7± ~ l7max =1=^01 atx = 7r/2 + n7r;a test charge comes to rest for 7q < 7^^^ twice per period, 
at phases that depend on u'q. We derive an explicit expression for the orbit of a particle in 
a LAEW with a triangular wave form in terms of elliptic integrals (Weierstrass functions). 
In the highly relativistic case, the exact solution for the orbit is well approximated by a set 
of light lines, with the velocity changing rapidly between ±c at phases x = 7r/2 + nvr. This 
leads to a triangular form for the orbit of a background particle, which oscillates about 
a fixed point. A test charge oscillates about a moving center, with the drift speed of the 
center of oscillation depending on /3q, but there is no simple analytic expression for the drift 
speed. 

Our motivation for this investigation was to explore the properties of linear acceleration 
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emission (LAE) in a LAEW, particularly the possibility of maser LAE in a LAEW as a radio 
emission mechanism, and high-energy LAE as an intermediate step in a pair cascade. The 
treatment of LAE inv olves some conceptua l and mathematical difficulties and is discussed 



separately in paper 2 (IMelrose fc Luo 



20091 ) ■ Here we write down the emissivity derived in 
paper 2, and use it to consider the absorption coefficient for LAE, with negative absorption 
implying maser emission. A mathematical difficulty is avoided by averaging the absorption 
coefficient over the small cone of emission. Maser emission is found to be possible for LAE 



in LAEW, as in earlier theories: 



Melrosd (119781 ) found maser LAE ir i a small ana plitude 



electrostatic oscillation with 7q ^ 7^3^^ in the present notation, and iRowd ( 119951 ) found 



maser 



(IRowe 



AE in a theory for LAE in a sinusoidal LAEW. We note that the latter theory 



19951 ) depends explicitly on a drift velocity, and in this sense is incompatible with 
the theory developed here and in paper 2. We discuss this difference and explain why there 
should be no explicit dependence on a drift velocity in LAE. 

Our estimate of the growth rate for maser LAE at low frequencies implies that it 
is too small for effective growth in a LAEW with an amplitude needed to trigger pair 
creation in an oscillating model for pulsars. For LAE in a LAEW to be viable as a radio 
emission mechanism, one requires a smaller amplitude LAEW, with 7max ~ 10- seems 
plausible that LAEWs with a wide range of amplitudes are present in different regions 
of the magnetosphere, and one needs to assume this to be the case for maser LAE to be 
effective. We conclude that maser LAE in a LAEW is a possible radio emission mechanism, 
but that the mere presence of a LAEW does not ensure effective maser LAE. We discuss 
the possible relevance of LAE to high-energy processes in pulsars in paper 2. 



We thank Matthew Verdon for helpful comments on the manuscript. 
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X X 

Fig. 3. — Trajectory of a test charge, on the left for 7q ^ 500, 7^^^^ ~ 1000, and on the right 
for Uq = 1.33 and for u^^x — ^ (lowest), 2.45 (middle) 1.25 (upper). As in Fig. [2], in the 
highly relativistic case the orbit is well approximated by a sequence of light lines. 
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A. Motion in a LAEW: covciriant treatment 

In this section we solve the equation of motion to find the orbit of a charge in a LAEW 
using a covariant formalism. 

This manuscript was prepared with the AAS lATgX macros v5.2. 
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A.l. Covariant formalism 

In the covariant formalism used here, greek indices run over 0-3, with signature 
+,—,—,—, and units with c = 1 are adopted. An electromagnetic field is described 
by the Maxwell tensor, F^'^, and for a LAEW with a phase which is an invariant, 
this is of the form F'^'^ = E{x)f^'^, where E{x) is the invariant electric field defined 
by 2E'^{x) = —F^^F^y. For a LAEW along the 3-axis, the tensor f^^ has components 
/"^ = —1, f^^ = 1, f^'^ = otherwise. The motion of the particle is described by its 
4-momentum = mu'^, where m is its rest mass and is its 4-velocity. Newton's equation 
of motion, in covariant form for the Lorentz force, is 

= qF^'^MuAr), (Al) 

where r is the proper time of the particle. The orbit may be written in the covariant form 

x'^ = 4 + A'^(r), «/^(r) = ^^J^, (A2) 

dr 

where x'^ = [t, x] is the 4-position and Xq contains the initial conditions. Our objective is 
to integrate (lAip once to find u^{t), and to integrate a second time to find X^{t). The 
LAEW is described by the wave form ([1]), with a wave 4- vector, = [17, K]. where 
Q is its frequency, and its wave vector, K, is assumed along the 3-axis. The invariant 
= K^K'' = ^2 _ |K|2 satisfied > for a superluminal LAEW. 

The one-dimensional motion is in the two-dimensional 0-3 subspace. The equation of 
motion then reduces to two invariant components, identifying by choosing two orthogonal 
4- vectors that span this subspace. We choose the 4- vectors K'^ and K'^ = f^^Ky, which 
are orthogonal, K'^K^ = 0. The invariants K^, Ku{t) = K^^u^It), Kou{t) = K^u^{t) are 
related by 

= -Kl = [Ku{t)]' - [Knu{r)]^ (A3) 
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The projection on the 4-velocity and the orbit onto these two 4-vectors gives 

Kujr) - KoujT) ET^ _ KXjr) K>^ - KoX{t) K^^, 

u [r) - , A [T) - — . (A4j 

respectively. An important step involves writing the phase of the LAEW in terms of these 
invariants. The phase is the invariant Kx, evaluated along the flA2p . so that the electric 
field in the equation of motion (lAip is that at the position of the particle. Ignoring a 
constant in the phase, it is identified as 

x(r) = KX(r), ^ = Ku{t). (A5) 

Using ( 1A5I) . the independent variable in (1A1I) is changed from r to x, and all quantities are 
regarded as functions of x- 

Projecting ( lAll) onto the 4-vectors K'^, gives 

dKu{T) .Kdu{t) dKDu{T) 

dx Ku[t) dx 

where we use ([I]), and where the dependence of x on r is implicit. 
Integrating the second of (]A6P gives 

Kdu{t) = KdUo + ujeF{x) , (AT) 

where KdUq is a constant of integration, and with 

F{x)=t / d^T{^). (A8) 

^0 



The solution of the first of (1A6P is 

Ku{t) = {K^ + [KDu{T)fY''', (A9) 
with only the positive square root allowed for K"^ > 0. 
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The function F{x) may be separated into a periodic part and a systematically 
increasing part, 



Fix) = Fix) + aoX, = ^ dxTix), (AlO) 

with Fix + 2tt) = Fix). The parameter is zero for any wave form that is symmetric, 
in the sense T(x + tt) = — T(x), and we restrict our discussion to such wave forms. The 
significance of ao 7^ is discussed in §5) 

The orbit is found by integrating dX^^ir) / dr = 'U^(r), which becomes 

dX^ir) _ Knuir) _ < + u^^jx) 

dx K^Kuir) ^ Kuix) ' ^ ' 

Integrating the first form in (lAllll gives 

where an invariant, dimensionless form of the displacement is given by 

71 \ ' KoUo + UJEFix') 







{K^+[Knu, + ujEFix')fYl 
where ( ]A7p and (IA9p are used. The proper time of the particle depends on phase through 



A. 2. Covariant form for the 4-current 



The source term for emission of radiation by any particl e is identified as the Fourier 



199l|). 



transform of the current density associated with the particle ( iMelrose fc McPhedran 
In a covariant formulation, this becomes the 4-current density 

J^'ik) = q [ druf^ir) e^^^M. (A15) 
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For motion in a LAEW, the only nonzero components of J^{k) are in the 0-3 plane. The 
formalism introduced in (IA4p allows one to write 



(A16) 



Using ( lA12p . the phase factor in ( 1A15I) becomes 



(A17) 



Using (jAlD, (1A15I) gives 



KJ{k) 


= (1 j dx 


1 




KDJ{k)_ 







(A18) 



with /3(x) = Kdu{x)/ Ku{x), and with now regarded as a function of x- 



The integral in (lAlSp may be evaluated after expanding in a Fourier series (IRowe 



19951 ). For a background particle, Z{x) is a periodic function and the expansion in Fourier 



series may be written 



1 



JZ{x)kKD/K' 



Us{k) 
Bs{k) 



(A19) 



with the Fourier coefficients given by 



Us{k) 
Bs{k) 



1 

2^ 



2-K 



dx 



1 



Jlsx+Zix)kKD/K^] 



(A20) 



Then flAlSP becomes 

KJ{k) 
KoJ{k)_ 

with the Fourier coefficient determined by 



KJs{k) 


= Q 


Us{k) 


KDJs{k)_ 




_Bs{k)_ 



KJs{k) 
KoJsik) 



27r5{s - kK/K"^). 



(A21) 



(A22) 
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For a test particle, one may write Z{x) = aiX + Z{x)i with a linear term, aix 
associated with the drift motion, and with Z{x) a periodic function. The Fourier series may 
be defined in terms of coefficients Us{k), Bs{k) defined by flA20p with Z{x) —>■ Z{x)- The 
result (]A22I) is modified by Us{k), Bs{k) Us{k), Bs{k) and s —> s — ai in the (5-function. 
After replacing the sum over s by an integral, performed over the 6 function, the result (!24l) 
is found to apply to both background and test particles. The separation of Z{x) into a drift 
term and an oscillatory term is irrelevant when treating LAE. 



